Infrared phonon activity in pristine graphite 
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We study experimentally and theoretically the Fano-shaped phonon peak at 1590 cm~^ (0.2 eV) in 
the in-plane optical conductivity of pristine graphite. We show that the anomalously large spectral 
weight and the Fano asymmetry of the peak can be qualitatively accounted for by a charged-phonon 
theory. A crucial role in this context is played by the particle-hole asymmetry of the electronic 
TT-bands. 



The discovery of graphene has triggered a tremendous 
amount of activity to study the electronic and struc- 
tural properties of single-layer and multi-layer graphitic 
compounds. The presence, in optical probes, of phonon 
anomalies related to the in-plane carbon displacements 
has in particular been useful both for the characteriza- 
tion and for the investigation of the fundamental prop- 
erties of graphene. While initially Raman spectroscopy 
was mainly employedft"— recent infrared experiments in 
bilayer and multilayer graphene have significantly con- 
tributed to the field, showing new remarkable phenom- 
ena, such as a strong modulation of the magnitude and of 
the lineshape of the optical phonon peaks as function of 
electrostatic gating^i^ the number of layerSfi^ and chem- 
ical doping.— 

The interest to graphene has also given a new impulse 
to the research on its 3D parent material, graphite. One 
of the long standing problems is the origin in graphite 
of the infrared phonon activity in the optical conductiv- 
ity at a; « 0.2 eV, which is associated with the in-plane 
Eu carbon displacements. Early measurements of the 
peak intensity, expressed in terms of the effective infrared 
charge Z, range from Z = 0.18 to Z = 0.41 electrons 
per atoms?i2iii^ qualitatively confirmed by more recent 
measurements fi^ii^ The possible origin of the infrared 
activity of the Eu mode in pristine graphite was theo- 
retically investigated in Ref. [l^ by means of a general- 
ized bond-charge model. The predicted phonon intensity 
was however significantly smaller than the experimental 
results. 

In this work we combine experiment and theory in 
order to gain new insights into this problem. Our in- 
frared measurements on highly ordered pyrolytic graphite 
(HOPG)^^ allowed us to characterize both the intensity 
and the Fano asymmetry of the phonon anomaly. We 
analyze these results within the context of the charged- 
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phonon theory,^'^ which has recently been shown to ac- 
count quantitatively for the phonon infrared intensity in 
few-layer graphenesJ^iii While in gated graphene the 
main factors that determine the infrared intensity are 
the induced doping and the internal electric field, here 
we show that in charge neutral graphite the phonon in- 
frared activity is essentially driven by the tight-binding 
terms breaking the particle-hole symmetry. 

Infrared reflectivity spectra were measured from 20 to 
6400 cm~^ at different temperatures from 10 K to 300 
K, using Fourier transform spectroscopy and in- situ gold 
coating as a reference. The spectra are shown in Fig. 
[T] The inset expands the region around 1600 cm^^, 
where the phonon peak is located. Optical conductiv- 
ity a{uj) was extracted by the Kramers-Kronig analysis, 
supported by spectroscopic ellipsometry between 6400 
and 36000 cm~^ as described in Ref. [l^. The spectral 
resolution was 5 cm~^. Fig. [^h- shows the real part of 
cf{uj) near the phonon peak, where a smooth baseline due 
to direct electronic interband transitions was subtracted. 
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FIG. 1: Infrared reflectivity spectra of highly ordered epitax- 
ial graphene at 300 K, 200 K and 10 K from RefUl. The inset 
shows the region near the phonon peak. 
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FIG. 2: (a) Real part of the optical conductivity Acri(aj) 
(thick black line) in the vicinity of the Eu phonon peak af- 
ter a smooth background is subtracted. Also shown is the 
fitting curve (red line) as obtained by Eq. ([l}. Panels (b)- 
(e) show the phonon parameters ujq, T, W, q extracted by 
the fit as a function of temperature (symbols) along with the 
theoretical calculations (dashed lines) not taking the phonon 
anharmonicity into account. 



The most striking feature is a pronounced peak asymme- 
try, signaling a Fano interference<^"— Accordingly, the 
spectra at every temperature are fitted very well with a 
standard Fano formulai^i^^ 



Aa'icj) 



2W +2qz-l 



(1) 



where z = 2(uj — loo)/T, and where ujq and F characterize 
the frequency and the linewidth of the phonon anomaly, 
while W and q represent the phonon spectral weight and 
the Fano asymmetry parameter, respectively. The tem- 
perature dependence of such phonon properties as ex- 
tracted from the fit procedure is shown in Fig. [^b)- 
(e). The phonon frequency decreases by about 2 cm^"'^ 
as the temperature increases from 10 K to 300 K. The 
spectral weight W increases from 950 to 1200 57~^cm~^, 
which corresponds to the variation of the infrared effec- 
tive charge e* = Ze from 0.28 to 0.31 e, in a good agree- 
ment with the previous experiments J^i^ii^ The mea- 
sured linewidth is close to 5 cm~^, which means that 
it is probably determined by the spectral resolution used 
and exceeds the actual value. 
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FIG. 3: (color online) (a) Atomic structure of graphite show- 
ing the relevant hopping terms ji. Atoms Bl and A2, con- 
nected by vertical 71, denoted by darker colors, contain also 
a local crystal field potential. 



The observation of a finite infrared phonon intensity in 
pristine undoped graphite is quite intriguing given that it 
is a monoatomic compound. No infrared intensity would 
be expected for instance if all the carbon atoms were 
equivalent. Moreover the evaluation of the static dipole 
in realistic graphite through a tight-binding calculation 
finds, like in bilayer graphene, a charge-disproportion 
between the inequivalent sites (and hence a static elec- 
tric dipole) three order of magnitude smaller than what 
observedi^ A suitable framework to account for the fi- 
nite infrared activity is thus the charged-phonon the- 
ory, which describes the borrowing of electronic effective 
charge by the optical mode from the electronic transi- 
tions to which it becomes coupledi^ Such a theory has 
been successfully applied to bilayer graphene to explain 
the strong dependence of the phonon intensity and Fano 
asymmetry as functions of gating. 

A key role in this context is played by the mixed re- 
sponse function x(a;),^^ which couples the current oper- 
ator with theelectron-phonon scattering operator related 
with the Eu lattice mode. The spectral weight is given by 
the relation W — 7r[Rex(wo)]^/'^o, while the asymmetry 
parameter is given hy q = — Rex(i^o)/Inix('^o)- 

The analysis in Ref. 17 was performed within a sim- 
ple tight-binding (TB) model, where only the nearest- 
neighbor in-plane and inter-plane hopping terms, 70 and 
71, were considered (see Fig. [3]). Such model can 
be generalized in a straightforward way for the bulk 
graphite by introducing the interplane dispersion, so that 
X(w) ^ X^°(a;) = /^^(dfc^/27r)x(fc^,w), where x{kz,uj) 
is the mixed response function for bilayer graphene with 
the replacement 71 — > 71(^2) = 271 cos{kz/2). Note how- 
ever that the description of the electronic structure of 
graphite in terms of this simple model would predict no 
infrared activity since the mixed response function would 
be zero. As we will show, the inclusion of high-order 
tight-binding terms, especially the ones that break the 
particle-hole symmetry, is fundamental to understand the 
origin of the phonon activity. 

In order to perform a quantitative analysis we consider 
the full tight-binding model defined in the 4x4 basis of 

(«Ik'^Ik'4k'^2k) where a|, 



the atomic orbitals ^I^j^ 



ik 
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and bl^^ represent the creation operators relative to the A 
and B sublattices in the layer i. The atomic basis and the 
different terms of the TB structure are shown in Fig. [3] 
The low energy bands of graphite around the K-H corner 
of the Brillouin zone are described by the Hamiltonian 
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where vr^ = fi{kx ± iky), 71.3^4 = 271 3 4 cos(fcz(i/2), 
72,5 = 272 + 272,5 cos{kzd), 1)3,4 = f73,4/7o- The param- 
agnetic current along the y-axis is thus evaluated as jy — 

I]k,<T *kJky*k where jky = -{e/h)dHw/dky. The cou- 
pled electron-phonon interaction operator can be written 
asV ^ g J2k <7 ^k^^k'/'y where g is the electron-phonon 
coupling constant as defined in Ref. [13, (t>y represents 
the dimensionless in-plane lattice displacement along the 
y direction. The precise form of V depends on the TB 
model considered. In a simple model, where only 70 and 
71 are present, Vk i n g raphite has the same expression 
as derived in Ref. |2l| for bilayer graphene. In a full 
TB model, however, Vk requires the knowledge of the 
derivative of 74 with respect to the lattice distortion u, 
in addition to the well known term djo/du. The matrix 
Vk thus reads 
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where F — cos(fczrf/2). The deformation potentials 
ao = d-ioldu = 4.4 eV/A ua = d'yA/du = 0.3 eV/A 
are extracted from LDA calculationsi^"— 

Once the current jy and the electron-phonon scattering 
operator V are defined, the phonon spectral properties 
can be evaluated in the framework of the charged-phonon 
theory following the derivation for bilayer grapheneJ^ 
We evaluate the mixed response function in the bare- 
bubble approximation, where we obtain the analytical 
expression: 



E 

k,n,m 



/(e«,k) - /(Cm.k) 



(4) 



where n, m are band indices, — (jka)"'"(Vk)™" and 
/(e) is the Fermi function. We introduced a phenomeno- 
logical damping term 77 to account for disorder effects. 
Finally, /3 — —igNyNg, where Ny = 2 and Ns — 2 stand 
for the valley and spin degeneracy, respectively. 

In contrast to the case of the simple 70-71 model, the 
inclusion of the higher-order TB terms makes xi'-^) de- 
pendent on the ultraviolet cutoff related to the finite 
bandwidth. To this aim we introduce in the linearized 
model a momentum cutoff k < k, where k is deter- 
mined by requiring the conservation of the Brillouin zone. 



TB model p-h symmetry W q 

full broken 666 -1.15 

partial: 72 broken 16 -0.53 
partial: 73 preserved 

partial: 74 broken 471 -0.94 

partial: 75 broken 50 -0.47 

partial: A broken 332 -0.74 

TABLE I: The spectral weight W and the asymmetry param- 
eter g at r = for the full TB model and for the partial TB 
models described in the text. W is expressed in units of Q,~^ 
cm~^. 



namely S'ceii = irNxk'^. We get thus k = 1.095 A"\ 
which is comparable with the distance F-K « 0.98 A~^ 
of the hexagonal Brillouin zone. 

To elucidate the role of the various TB parameters on 
the infrared phonon peaks, we consider for the moment 
the full set of TB parameters from Ref. [2^. The phonon 
spectral properties obtained for this model retaining all 
TB parameters, using rj = 5 meV and T=0 K, are re- 
ported in Table HI As we can see, the inclusion of the 
higher-order terms is crucial to give rise to a finite in- 
frared activity. 

In order to investigate the role of each TB parame- 
ter on the infrared activity, we consider a set of toy TB 
model, based on our present parameter set where, keep- 
ing the terms 70 and 71 fixed, we "switch on" , one by 
one, the remaining TB parameters 72, 73, 74, 75 and the 
crystal field A. The chemical potential is evaluated self- 
consistently in each model to preserve the charge neutral- 
ity. The results are summarized in Table U where one can 
clearly see how the onset of infrared activity is directly 
related to the breaking of the particle-hole symmetry. 
Within this scenario, the detailed optical properties re- 
sult from a complex balance between all TB parameters. 

To estimate the sensitivity of the predictions of the 
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FIG. 4: Phase diagram of the phonon strength W vs. the 
Fano parameter q. The empty square represents our experi- 
mental values, filled circles represent data obtained at T 10 K 
by using different TB models available in literature.—^ The 
size of symbols for the theoretical data is proportional to the 
damping term 77 = 5, . . . , 30 meV. 
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charged-phonon model on the particular choice of the 
TB model, we calculated the phonon strength W and the 
Fano factor q obtained from different TB parameter sets 
proposed in literature for graphite. ^^^-^ The results are 
summarized in Fig. |4j where we also show the experimen- 
tal values of W and g at 10 K. The large spread of the the- 
oretical points in the W-q diagram reveals that relatively 
small changes in the underlying electronic structure af- 
fect significantly the phonon properties. Given the high 
sensitivity of the results of the charged-phonon model on 
the electronic details and some uncertainty related to the 
choice of the momentum cutoff fc, the comparison with 
the experimental finding seems to be reasonably good. 
Importantly, the model predicts correctly the large value 
of the experimental spectral intensity W and therefore 
the effective infrared charge. This is a significant im- 
provement as compared to a previous calculation based 
on a generalized bond-charge model,-^ which predicts the 
effective charge of this mode to be an order of magni- 
tude smaller. The charge-phonon theory also accounts 
for (and even overestimates) the strong Fano lineshape 
asymmetry. It should be noted, however, that the pre- 
cise value of the Fano parameter q is also quite sensitive 
to the damping magnitude, which is sample dependent 
and should be regarded as an adjustable parameter in 
the present model. 

It is interesting to compare the temperature depen- 
dence of the phonon parameters shown in Fig. [2] with 
the theoretical calculations (dashed lines). The frequency 
shift and the linewidth of the phonon peak here are 
evaluated via the phonon self energy 11 (cj) : Aluq — 
h-^Re[n{ujo)],T = -h~Hm[Il{ujo)], while the spectral 
weight and the asymmetry parameter are obtained us- 
ing the given above formulas. We can notice that for all 



parameters the theory strongly underestimates the tem- 
perature variation as compared to the experiment. We 
believe that the strong temperature dependence of the 
experimental features is a signature of anharmonic ef- 
fects induced by the phonon-phonon scattering. A care- 
ful study of such anharmonic effects within the density 
functional theory context was provided in Ref. 0, where 
these affects were found to be mostly responsible for 
the thermal variation of the phonon frequency and the 
linewidth. It is likely that the phonon anharmonicity is 
also responsible for a sizeable temperature dependence of 
the spectral weight and the Fano asymmetry. The the- 
oretical treatment of anharmonicity within the charged- 
phonon theory is an interesting problem which is beyond 
the scope of the present work. 

In conclusion, we have investigated the origin of the in- 
frared phonon activity of the Eu mode at w w 1590 cm~^ 
in pristine graphite. We experimentally found that the 
effective infrared charge is close to 0.3 electrons, which is 
an exceptionally large value for an elemental compound. 
A significant Fano asymmetry is also present {q w -2). 
We have shown that both these observations can be ex- 
plained within the framework of the charged-phonon the- 
ory, where a crucial role is played by the broken particle- 
hole symmetry associated with the higher-order tight- 
binding terms. The approach used in this work can be 
extended to other related materials, such as carbon nan- 
otubes and intercalated graphite. 
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